Chaotic quantum ratchets and filters with cold atoms in optical lattices: properties of 

Floquet states 
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Recently, cesium atoms in optical lattices subjected to cycles of unequally-spaced pulses have 
been found to show interesting behavior: they represent the first experimental demonstration of a 
Hamiltonian ratchet mechanism, and they show strong variability of the Dynamical Localization 
lengths as a function of initial momentum. The behavior differs qualitatively from corresponding 
atomic systems pulsed with equal periods, which are a textbook implementation of a well-studied 
quantum chaos paradigm, the quantum 5-kicked particle (5-QKP) . We investigate here the properties 
of the corresponding eigenstates (Floquet states) in the parameter regime of the new experiments 
and compare them with those of the eigenstates of the <5-QKP at similar kicking strengths. We 
show that, with the properties of the Floquet states, we can shed light on the form of the observed 
ratchet current as well as variations in the Dynamical Localization length. 

PACS numbers: 32.80.Pj, 05.45. Mt, 05.60.-k 



I. INTRODUCTION 

Periodically- kicked quantum systems, such as the 5- 
kicked particle (<5-KP), have long played a central role 
in studies of quantum chaos and the correspondence 
between quantum behavior and the underlying classi- 
cal dynamics 0, 0. Advances in the manipulation of 
cold atoms have permitted the experimental realization 
of these systems in pulsed optical lattices Experi- 
ments with sodium and cesium atoms have demonstrated 
the phenomenon of "Dynamical Localization" (DL) - the 
quantum suppression of classical chaotic diffusion - and 
established the suitability of these systems as an arena 
for the study of effects arising from quantum chaos. 

When treating conservative quantum systems it is fre- 
quently useful to study the system's energy-eigenstates, 
and for periodically driven systems the appropriate gen- 
eralization of these states is given by the Floquet states. 
This approach has provided extensive insight into the 
properties of the standard QKP, and has shown, for ex- 
ample, that DL arises directly from the exponential lo- 
calization of the system's Floquet states 0- Observed 
momentum oscillations associated with chaos-assisted 
tunneling, in experiments using periodically-driven cold 
atoms Q and BECs [f| have also been analysed with 
Floquet theory; it was found that the oscillation period 
is determined by the splittings of the Floquet phases of 
a pair of symmetry-related eigenstates. The statistics of 
QKP Floquet quasi-energy spectrum have been studied 
extensively see e.g. |(| and compared with the predictions 
of Random Matrix Theory. Notably, though, the <5-QKP 
has Poissonian short-range statistics (which are typically 
associated with integrable dynamics) even for very large 
values of K, where the dynamics is fully chaotic. This 
has been shown to be a further consequence of DL |(| . 

However, a series of recent theoretical 0, 01 an d ex- 
perimental 0, 0, studies of cold atom systems sub- 
jected to repeating cycles of unequally spaced kicks re- 
vealed dynamics rather different from that found in the 



corresponding standard QKP systems. Two types of un- 
equally (5-kicked systems were investigated. The first 
consists of a <5-KP with small perturbations made 
to the kick-period. We term it the perturbed-period KP. 
In the second system, the system is periodically sub- 
jected to pairs of closely-spaced kicks. This is referred to 
as the double S-KP or 26-KP. 

In a theoretical study, the perturbed-period KP was 
found to yield a quantum momentum current even in 
the chaotic regime Q, 0- This was unexpected in a 
chaotic Hamiltonian system, since to date only mixed 
phase-space ratchet mechanisms had been investigated 
[l2lll3l| . A simple definition of a ratchet is a spatially pe- 
riodic device which produces a current without net bias 
(ie the time and space- averaged forces are zero). Most of 
the extensive ratchet literature deals with dissipative or 
Brownian ratchets 01 an d comparatively little theoret- 
ical work has been undertaken on Hamiltonian ratchets, 
which are dissipation and noise- free. In |gL IXdl| . a momen- 
tum distribution with a non-zero average (constant in 
time) was obtained experimentally from an atomic cloud 
with initial zero average momentum. We are unaware 
of any other experimental studies of Hamiltonian quan- 
tum ratchets: all implementations to date have been of 
dissipative/Brownian ratchets. Hence the results from 
|10| and also reproduced here, represent the only imple- 
mentation of a type of Hamiltonian quantum ratchet, 
whether chaotic or mixed phase-space. In Ref. || it 
was also proposed that the chaotic diffusive properties 
of the perturbed-period KP could be exploited to filter 
cold atoms, according to their momenta, by controlling 
the Dynamical Localization. 

For the second system, the 25-KP, a theoretical and 
experimental study revealed that the diffusion is 
dominated by long-ranged correlations which control es- 
cape from well-defined momentum trapping regions. This 
combination of strong chaotic diffusion and long-ranged 
correlations is novel: strong chaos is generally associated 
with rapidly decaying correlations. 
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It is clear that Floquet theory is central to the analy- 
sis of chaotic, time-periodic quantum systems. The need 
to understand further the chaotic Hamiltonian ratchet 
as well as the 25-KP motivated this study of the Floquet 
states of these systems. The paper is organized as follows. 
In the next section we review the well-known <5-KP, then 
introduce the perturbed 6-KP and the double S-KP sys- 
tems. In Section III we give a brief review of the Floquet 
approach. In Section IV we compare the results with re- 
cently obtained experiments on these systems. Finally in 
Section V we give our conclusions. 

II. INTRODUCTION TO <5-KICKED SYSTEMS 

The Hamiltonian for the usual <5-KP can be written as 

2 

H =y + Ksmx^28(t-nT) (2.1) 

n 

where K is the kick strength and T is the time interval 
between successive kicks. Consider its effect on an ensem- 
ble of particles with a gaussian momentum distribution 
centered on p , N(p) = exp[— (p — p ) 2 /Ap 2 }. The clas- 
sical dynamics depends only on the parameter K, and 
for values larger than K ~ 1, the chaotic diffusion is not 
bounded by classical barriers. In this regime the ensem- 
ble will diffuse in momentum space, its average energy 
growing linearly with time as (p 2 ) = Dt, where, to low- 
est order, the diffusion rate is given by Dq = K 2 /2. The 
distribution will thus remain gaussian, although its width 
will increase with time as Ap(t) = \J Dt. 

In contrast, the quantum system only follows this be- 
havior up to a timescale t* ~ D/h 2 |15|. after which 
diffusion in momentum space is suppressed - dynami- 
cal localization (DL). Such a system will asymptotically 
evolve towards a characteristic exponential momentum 
distribution, N(p) ~ exp [—\p — po\/Apq], with constant 
width Apq ~ \jDt* ~ D/h, which thus acts as an exper- 
imental fingerprint for DL [3( . As DL is a wave-coherent 
effect, the quantum system must preserve coherence over 
at least the timescale t* for this effect to be observable. 
Indeed, it has been verified experimentally ^(| that the 
DL profile does not survive the presence of noise or dissi- 
pation, and that with decoherence a more gaussian profile 
for N(p) will be produced. 

In both classical and quantum cases the behavior of the 
standard i5-KP is essentially independent of po since, even 
for modest values of K, the effects of small fluctuations 
in the structure of phase-space are on negligible scales 
relative to ApQ. Even if there are small stable islands, 
they are of size Ap ~ 1 so have little effect on the general 
form of N(p), since typically ApQ ^> 1. 

The classical dynamics of the <5-KP is obtained by 
iterating the well-known "Standard Map". For the 
perturbed-period and double 8- kicked systems, on the 
other hand, the dynamics is given by a 2-kick map: 

Pj = Pj-i ~ v '( x j-i) 



Xj — Xj-i + PjTi 
Pj+i = Pj ~ V'(xj) 
Xj+i = xj +p j+ xT 2 . 

Clearly, setting T\ = T 2 and V'(x) = Kcosx, we recover 
the Standard Map. For the perturbed-period KP, the 
lengths of the two kicking periods are T\ = 1 + e and 
T 2 = 1 — e, where e< 1. The perturbation thus consists 
of slightly altering the kicking period about its mean. For 
the double S-KP we take T x = 2 - e, T 2 = e, although 
we shall also show the effect of interchanging the two 
kick periods. It should be noted that these systems are 
time-periodic, with period T to t = T\ + T 2 , and are thus 
quite distinct from the recent interesting study of two in- 
dependent kicking sequences, which can be non-periodic 
and hence non-localizing |17|. 

As in the standard map, we consider a sinusoidal po- 
tential V(x) = if sin 2;. However, to obtain a ratchet 
current in the case of the perturbed-period KP, we need 
to break the spatio-temporal symmetries, and so we 
add an additional "rocking" linear potential of strength 
A. In this case the form of the potential is V(x) = 
— [Ksmx + Ax{— 1) J ], where j is the kick number. In 
experimental implementations of this system, the rock- 
ing linear term was obtained by means of an accelerated 
lattice Q. 

We first consider in general terms how the introduction 
of the second timescale e modifies the classical behavior 
of the standard map. If we neglect all correlations, the 
standard map has a constant momentum diffusion rate, 
Do ~ K 2 /2 - this is what one would expect if the mo- 
menta at consecutive kicks are uncorrelated and so evolve 
as a random walk. However, unless K is exceedingly 
large, the time-evolution of the standard map will con- 
tain some short-range (2-kick and 3-kick) correlations. 
Including these corrections yields a modified diffusion 
constant D = ^-[1 - 2J 2 {K) - (Ji(K)) 2 . . .}. Of par- 
ticular interest is the J 2 {K)K 2 term, representing corre- 
lations (V'(xj)V' (xj+2)) between nearest-but-one kicks 
(the 2-kick correlation). 

For the modified systems it is also possible to obtain 
analytically the important correlations @, EJ- F° r m_ 
stance, for the perturbed-period system the corrected dif- 
fusion is D ~ ^[l~2J 2 {K)cos{2p e^A)-{J l {K)) 2 . . .]. 
We see that we have a modified 2-kick correlation which 
oscillates with the initial momentum pq. This effect 
is clearly most significant for values of K such that 
2J 2 (K) - 1. 

The key point is that perturbing the kick spacings T 
by a small amount can result in large scale (relative to 
ApQ ~ D/h) variations in the classical momentum dif- 
fusion, and that these are present even in fully chaotic 
regimes (we take this to mean the absence of visible sta- 
ble structures on the Poincare surface of section). For 
the analysis of experiments, one must now consider a 
local diffusion rate D(po), which depends on the initial 
relative momentum between the atoms and the optical 
lattice. 
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In |9(, the perturbed-period system was implemented 
experimentally with a cloud of cesium atoms for the case 
A = 0. It was verified that the energy absorbed by 
the cloud after Dynamical Localization ((p — po) 2 ) oc 
cos(2poe) as expected. However, A — corresponds to a 
symmetric potential. The case A — n/2 in the perturbed- 
period system is particularly interesting since then the 
momentum diffusion is asymmetric about po = 0. This 
implies that atoms with positive momenta will absorb ki- 
netic energy at different rates from those with momenta 
of the same magnitude but moving in the opposite direc- 
tion. This asymmetric momentum diffusion represents a 
type of fully chaotic momentum ratchet: in other words, 
roughly equal numbers of particles will diffuse to the left 
or to the right, but those diffusing to the right, on aver- 
age, move faster. 

In a further experimental study of the perturbed- 
period KP ^3]> the A — 7r/2 potential was implemented 
by means of an accelerated lattice. It was found that 
an atomic cloud prepared initially with a gaussian mo- 
mentum distribution centered on po — evolved into a 
distribution with non-zero, but constant, (p) which per- 
sisted even beyond the break-time t*, as expected from 
the theory 0. This type of chaotic directed motion was 
first identified in a slightly different system: a kicked 
asymmetric double- well potential Q- However, the latter 
potential gives rise to a rather more complicated diffusive 
behavior, and also has proved much harder to implement 
experimentally. For these reasons here we do not con- 
sider the case of the asymmetric double- well ratchet, but 
note that our Floquet analysis of the perturbed-period 
KP can be carried-over to the system investigated in • 

The second system we consider explicitly in this work, 
the 25-KP, has diffusive behavior which is qualitatively 
different to both the standard map and perturbed-period 
KP. While for these other kicked systems we can analyze 
the diffusion as an uncorrelated term, K 2 /2, corrected 
by short-ranged correlations (typically only 2 or 3-kick 
correlations for K ~ 3), for the double i5-KP, we find 
that the diffusion at long times is dominated by families 
[TTj l of long-ranged "global" correlations ( "global" in the 
sense that they correlate all kicks up to the time under 
consideration). At short times, the diffusion is dominated 
by a 1-kick correlation not present in other kick systems. 
At longer times, the global diffusion terms, though weak, 
accumulate and eventually become dominant. 

The method of correlations provides a generic and ac- 
curate way of interpreting experimental data for this sys- 
tem 0| • There is also a simple physical picture. For 
particles subjected to kicks of form i'Tsina;, consecu- 
tive kicks will be out of phase and will hence cancel if 

poe ~ (2n + l)ir where n = 0, 1, 2 In other words, 

an impulse V'(xj) = Ksmxj will be immediately fol- 
lowed by another which cancels it, since V'(xj+i) = 
K siriXj+i ~ K sm(ir + Xj). This cancellation means that 
particles become trapped at these momenta. In contrast, 
particles for which p e ~ 2n7r, will experience enhanced 
diffusion. 



It was shown in that the new types of global fami- 
lies of correlations control the escape from, and through, 
these "trapping regions" . An unexpected feature of the 
classical calculations (and also seen in experiment) was 
the observation that particles initially prepared in the 
trapping regions will eventually gain more energy than 
those initially prepared in regions of enhanced diffusion, 
after a timescale t 3> l/(JtTe) 2 [lj. 



III. QUANTUM DYNAMICS AND FLOQUET 
STATES 



If a Hamiltonian has a T-periodic time dependence, 
H(t + nT) = H(t), then the Floquet theorem implies 
that solutions to the time-dependent Schrodinger equa- 
tion can be written in the form 



ip(t) = exp [-iet/K\ (j>{t) 



(3.1) 



where 4>(t) is a T-periodic function called a Floquet state, 
and e is a quantity with dimensions of energy, termed 
a quasienergy. This type of relation is familiar in the 
context of solid-state physics, where a Hamiltonian's in- 
variance under discrete shifts of spatial position (typi- 
cally arising from a lattice structure) allow solutions to 
be written analogously in terms of Bloch states and quasi- 
momenta. The Floquet states provide a complete basis, 
and thus the time-evolution of a general quantum state 
under periodic driving can be expressed as 



= 2_j C « eX P [~ ie nt/K\ (f> n (t), 



(3.2) 



where {c„} are time-independent expansion coefficients. 
It is clear from this expression that the Floquet states and 
quasienergies play a similar role for periodically-driven 
systems to that of energy eigenstates and eigenvalues in 
the time-independent case. 

The time-evolution operator U{t2,t\) may be used to 
evolve a quantum state from time t = t\ to a time t = ti . 
For a time-periodic system, the single-period propagator 
U(T, 0) allows a quantum system to be evolved "strobo- 
scopically" at intervals separated by the period T with 
great efficiency, by defining the quantum map f (nT) — 
[U(T,0)] n *(0). In terms of Floquet states, it is straight- 
forward to show that the time-evolution operator is given 
by 



U(t 2 ,h 



Texp 



H(t')dt' 



(3.3) 



J2^ £nt/h \Ut2))(Mh)\, (3.4) 



where T is the time-ordering operator, and thus it can 
be seen that the quasienergies and Floquet states can be 
conveniently obtained by simply diagonalising the one- 
period propagator. The eigenvectors of this operator are 
the Floquet states, while its eigenvalues are related to 
the quasienergies via A„ = exp [— ie n T /%]. 
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Obtaining an explicit form for this propagator is nor- 
mally a complicated procedure, as in general the driving 
field does not commute with the static Hamiltonian. For 
the case of 5-kicking, however, the problem is simplified 
considerably which allows an analytic form for the prop- 
agator to be written. For the QKP fEa 12. 1|) the propa- 
gator is given by 

U(T, 0) = exp [-iTp 2 /2h] exp [-i(K/K) sinx] . (3.5) 

Using a basis of plane-wave states, the matrix-elements 
(rn+q\U(T, 0)\n+q') of this operator can easily be shown 
to be 

U m , n (T, 0) = exp [-iT(m + q) 2 /2h] J m _ n (K/h)5(q-q'), 

(3.6) 

where q is the quasimomentum (p = (m + q)h) and J n is 
the n-th Bessel function of the first kind. For practical 
purposes it is useful to note that |J m -„(a;)| decreases 
extremely rapidly with increasing \m — n\, thus giving U 
an effectively banded-structure. 

The single-period propagators for the unequally-kicked 
systems can now be expressed as the product of two ma- 
trices of this form, U(T tot ,0) = U{T X + T 2 ,T 1 )U(T 1 ,0), 
where T\ and are the two kick-periods. It should be 
noted that although the single-kick propagators do not 
conserve quasimomentum, their product does. As a con- 
sequence the evolution of an ensemble of non-interacting 
particles which can be modelled by the evolution of a su- 
perposition of states with different quasimomenta, repre- 
sents a computationally efficient procedure: we can con- 
sider each quasi-momentum component independently. 
In our study of Floquet states in effect this means we 
can diagonalise the smaller matrix given bv 13.61 which is 
block-diagonal in q. 

Having obtained the Floquet states by diagonalising 
(|3.6|) . it is useful to analyse their structure, particularly 
their spread in momentum space. We do this by evalu- 
ating a localization length, L of each Floquet state at 
t = 0. We note that in this case we cannot assume 
that the Floquet states have the usual exponential mo- 
mentum distribution N(p) ~ exp [— \p — po\/L], of the 
usual 5-KP. Hence we take L to be simply the root mean 

square deviation from the mean, L = \J p 2 n — p^. In 
this expression p n is the mean momentum of the n- 
th Floquet state at t = 0, p n = (<j) n (0)\p\(f) n (0)) , and 
p 2 n = {(j) n (0)\p 2 \(j) n (0)) . If the Floquet states do not have 
a strong time-dependence, such as for the standard, or 
even the perturbed-period KP, this is adequate to quan- 
tify the degree of spreading in momentum space. We 
shall see, however, that for the case of the 28 kicked sys- 
tem it is not sufficient to measure the localization at a 
single time, due to the extremely strong time-dependence 
of the Floquet states. 



IV. RESULTS 
A. Perturbed-period KP: a chaotic ratchet 

FigQ] shows a plot of the experimental ratchet current 
obtained in [Toj |. A series of momentum distributions 
N{p — Pa) as a function of po for a cloud of cold cesium 
atoms in an optical lattice pulsed with unequal periods. 
The momentum distribution is essentially unchanged af- 
ter about 60 kicks; the plotted values correspond to about 
T = 200 kicks, hence well after Dynamical Localization. 
Full details are given in [jj 0] , but by employing an ac- 
celerated lattice, the experiment simulated an effective 
rocking potential with A ~ 7r/2. The first moment of 
each localized distribution I = (p — po) was then calcu- 
lated and plotted as a function of po- In particular, a 
distribution centered at p = initially, and with zero 
initial momentum current (p) = at t = 0, yielded a 
finite and constant momentum current (p) ~ 4 at long 
times. For non-zero initial momenta po ^ 0, an oscilla- 
tion / cx cos(2poe), was observed and is seen in FigQ 
This may be qualitatively understood from the form of 
the classical two-kick momentum-diffusion correction in- 
troduced in Sec. II, C*2 = — K 2 cos(2poe — ir/2). If we con- 
sider a very small momentum displacement 8p = p — p , 
differential absorption of energy for particles moving to 
the left 5p < or right 5p > is proportional to the gra- 
dient cx K 2 cos (2poe). An accurate analytical form 
for the classical current was derived in 

Of course, it follows that a finite and persistent con- 
stant momentum current is also obtained classically. It 
was found in Q that asymmetric diffusion persists only 
on a timescale t ~ l/(Ke) 2 and for this (unbounded) 
chaotic system, the acquired momentum asymmetry is 
never lost. For a bounded ("compact phase-space") sys- 
tem, such asymmetries would vanish on a long time-scale, 
since the distribution of a fully chaotic system would 
eventually become uniform. For this reason, until re- 
cently, it was argued that a fully chaotic system could 
not generate directed motion. So, although as shown in 
, the fully chaotic classical system can keep a constant 
current for long times, practical implementation is ulti- 
mately less interesting since the average kinetic energy of 
the ensemble grows linearly with time and without limit. 
Hence, this type of chaotic ratchet is of most interest as 
a quantum rather than a classical ratchet since in the 
quantum case DL halts the diffusion and "freezes-in" the 
asymmetry, without the need for classical barriers like 
tori. 

In Fig[5] we reproduce two experimental momentum 
distributions for K ~ 3 obtained with cesium atoms in 
Rcf. 9], for A = ±7r/2. We clearly see that the origin of 
the non-zero momentum current is the asymmetric mo- 
mentum distribution. 

As expected, Fig|2 shows that changing the sign of A 
reverses the asymmetry. At this stage it may be unclear 
to the reader what the significance of altering the sign of 
A in the experiment might be, since after all, the rocking 
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□ Experiment: K=2.6 E =1/16 A=7i/2 k 1 
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FIG. 1: Experimental values of the momentum current I = 
{p — po), for the perturbed-period KP, obtained with cold 
cesium atoms in a pulsed optical lattice for K ~ 3, e — 1/16. 
The solid line is a best-fit to the data, showing that the current 
oscillates sinusoidally as / <x cos2po£ for A — tv/2. 




FIG. 2: Experimental momentum distributions N(p) for the 
perturbed-period KP obtained with cold cesium atoms in a 
pulsed optical lattice for K ~ 3, e = 1/16 . The distributions 
have localized, and hence remain essentially constant with 
time. The results show clearly that the origin of the net non- 
zero value of (p) obtained at long times is in the asymmetry 
of the DL profiles. As expected, the asymmetry is reversed by 
changing the sign of A, the amplitude of the rocking potential. 



potential involves alternating impulses if since ± A. In 
fact the distinction (as may be ascertained from the form 
of the classical diffusion) is between the case where an 
impulse K sin x + A precedes free evolution for a time 
interval T\ = 1 + e (obviously followed by an impulse 
if sin a; — A and interval T2 = 1 — e) and the separate 
experimental case where an impulse K sin x — A precedes 
free evolution for a time-interval T± = 1 + e and so forth 



(which corresponds to a reversed current) . 

Note that the experimental range of K ~ 2.6 — 3.4 
does correspond to a classical surface of section with 
some islands. However, we note that classical quanti- 
ties such as the average energy are very accurately given 
by diffusion rates (with 2 and 3-kick corrections). The 
essential mechanism is asymmetric chaotic diffusion: sim- 
ilar behavior was found at larger K in ,7| in regimes 
where there are no visible classical islands (but for which 
experiments are not available); hence, in the analysis 
of this type of ratchet, the presence (or otherwise) of 
small stable islands is immaterial. What is important, 
though, is that since the asymmetric diffusion term is 
2J 2 (if)cos(2p e - A), we need J 2 {K) ^ 0. Thus the 
much-studied (for the standard map) parameter value 
K = 5 does not produce asymmetry, since 0/2(5) — 0. 
Values of K ~ 2.5 — 3.5, % = 1/4— 1, on the other hand, 
turned out to be experimentally convenient and produced 
the strongest asymmetries. 

11)0 
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FIG. 3: The graph shows the localization lengths L of Flo- 
quet states as a function of average momentum p (i.e. (p)) for 
K — 3.4, h—1. Results are shown (a) for the standard QKP 
case, i.e. A = 0.01 (a non-zero A was used to break spatial 
symmetry), e = 0, (b) for e = 0.01, A = 0.01 and (c) e = 0.01, 
A — tt/2. The graph shows that for the standard kicked rotor 
the L are distributed within a narrow range in comparison 
with other two below. For the rocking case, L oscillates with 
p as expected from the 2-kick correction 2J2(K) cos(2poe — A); 
the oscillations of the two lower graphs are shifted relative to 
each other by a phase n/2. The density of eigenstates corre- 
sponding to average momentum range is roughly the same in 
all three cases. 

We now examine the form of the underlying Floquet 
states. In Fig|3]we compare the localization lengths for 
the standard QKP, with those of the perturbed-period 
KP for K — 3.4, e = 0.01. The difference is quite striking; 
while the standard QKP eigenstates are quite uniform 
across all regions of phase-space, the perturbed-period 
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localization lengths oscillate sinusoidally with p, with a 
period of 7r/e. Introducing the additional rocking poten- 
tial with the accelerated lattice (A = it/ 2) clearly leads 
to a 7r/2 shift in the oscillations. Inspecting Fig[3{c) for 
p ~ 0, we see that for positive momenta the localiza- 
tion lengths are increasing, while for negative momenta, 
the localization lengths decrease. Note the nearly regu- 
lar row of states for the standard QKP case with L ~ 1 . 
These correspond to states localized on a series of stable 
islands separated by 2tt, due to the momentum periodic- 
ity of phase-space in that case. 

We have chosen a parameter range for which L -C ir/e: 
that is, the localization length of each state is much 
smaller than the oscillation in p. Hence individual Flo- 
quet states really do sample "local" diffusion rates. We 
found that if we move towards a regime where L ~ 7r/e, 
the conclusions remain valid, but the amplitude of the 
oscillations is considerably damped. Similarly, if the sign 
of J%{K) changes, so does the sign of the sinusoidal os- 
cillation. 

We now consider the shape of the Floquet states in 
detail. In FigH we show the momentum distributions 
N(p) = \(f> n (p)\ 2 for Floquet states of the standard QKP. 
The distributions (with N(p) on a logarithmic scale) all 
show the well-known triangular form Q - the hallmark 
of dynamical localization. It may be clearly seen that the 
localization lengths vary little from state to state. 
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FIG. 4: Floquet states for the standard QKP, for K = 3.4, 
h — 1. As expected, all the states are exponentially local- 
ized, giving the characteristic triangular shape of N(p) when 
plotted on a logarithmic scale. They all have approximately 
similar localization lengths. 
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In FigO by contrast, the localization lengths of the 
Floquet states of the perturbed period <5-KP display a 
strong dependence on the mean momentum of the states. 
In addition, the figure shows that states localized close 
to p = are markedly asymmetric. The states are con- 
siderably extended towards positive momentum, but are 
strongly localized towards negative p. This behavior 
neatly accounts for the form of the experimental momen- 
tum distribution shown in Fig.l, which for A = n/2 were 
also more extended towards positive p. The states local- 
ized near p ~ 7r/4e and 7r/4e correspond to, respectively, 
minima and maxima of the classical diffusion. They are 
roughly symmetrical (typically) but vary by up to a fac- 
tor of ~ 40 in L. In Ref.@ it was proposed that the 
observed variation in the energy absorption rates be- 
tween atoms prepared with an initial drift momentum 
Po = — 7r/4e (which absorb very little energy) and those 
with po — ir/4e might be exploited to filter traffic of 
atoms through an optical lattice. The form of the un- 
derlying Floquet states explains this differential rate of 
energy absorption. 

Subsequently, it was found experimentally that the 
double <5-KP in fact shows much more pronounced dif- 
ferential absorption rates, without requiring the applica- 
tion of a rocking field A. We next report a study of the 
Floquet states of this system. 



FIG. 5: Typical form of Floquet states for the perturbed- 
period KP, K = 3.4, e = 0.01, A = tt/2 and ft=l. Here we 
plot N(p) = |0 n (p)| 2 as a function of p. (a) states with p ~ 
— 78. This corresponds to a minimum of the 2-kick correction 
cos2poe — tt/2. The states are narrow, but in general, roughly 
symmetric, (b) states with p ~ 0. The typical state here is 
asymmetric (c) states with p ~ +78. This corresponds to a 
maximum of the 2-kick correction. States here are generally 
symmetrical, but broad and flat-topped. 



B. Double (5-KP 

A study of the experimental and classical behavior of 
the double (5-KP was carried out in 11]. The classical 
dynamics is very different to that of the perturbed pe- 
riod KP. At very short times, the chaotic diffusion com- 
prises an uncorrelated diffusion term K 2 /2 and one dom- 
inant 1-kick correction. It was found in pd| that one can 
approximate the growth in the mean energy with time 
t, by the simple expression (p 2 ) ~ K 2 t[l — cospo £ ]- In 
Fig|SJa) we show experimental results for cesium atoms 
which localized in this regime. The experiment measured 
the energy of a series of clouds of ~ 10 6 atoms moving 
through the pulsed optical lattice with varying average 
drift momenta pq. For Fig|Sfa), the simple expression 
given above gives an excellent fit to the experiment, if 
we take t ~ t*, where t* is the break time. This regime 



7 



600 



500 
A 400^ (a) 
0,300 
v 200 

100 




£ = 0.045 




FIG. 6: Experimental results for double S-KP realization with 
cesium atoms (see [Tl| for details). Each data point (star) 
shows the energy absorbed (after 100 kicks, K = 3.3, fi — 1) 
by a cloud of atoms with average momentum p — po (rela- 
tive to the optical lattice) at initial time, t = 0. With in- 
creasing e, we see the minima (maxima) in the energy flip 
into maxima (minima) as a long-ranged family of classical 
correlations gradually overtakes the 1-kick classical correla- 
tion. The dashed lines represent a classical simulation using 
100,000 particles, all with momenta po at t — 0, and K within 
the range 3.3 ± 10%. (a) t* < ti ~ l/{Kef. Regime where a 
one-kick correlation is the dominant correction to the classi- 
cal diffusion. Here, atoms prepared near the trapping regions 
(poe ~ (2n + l)ir) remain trapped. Results follow closely the 
formula (p 2 ) ~ K 2 t(l - cosp e). (b)i* ~ l/(Ke) 2 . Regime 
showing the inverted peaks of the Poisson correlation terms 
analyzed in [Tl|. which determine the momentum trapping 
very close to the resonant condition (poe = (2n + l)ir). (c) 
t* > l/(Ke) 2 . Regime dominated by correlation family Cgi, 
but sharp inverted peaks due to the Poisson correlations are 
still visible. 



corresponds to t* <C l/(Ke) 2 . 

However, a more detailed study of the classical corre- 
lations showed that for later times, a new type of cor- 
rection appeared. Families of long-ranged, or "global" , 
correlations which coupled all kicks appeared. These cor- 
rections are individually very weak, but accumulate to 
eventually dominate the diffusive process. One family 
(termed the "Poisson family" in 1J|) was shown to lead 
to well-localized, inverted peaks in the energy absorp- 
tion at values of po ~ (2n + l)%/e, where n = 0, 1, 2... . 
These values of po correspond to trapping regions in 
phase-space (at low values of K, structures correspond- 
ing to islands and broken phase-space barriers are evi- 
dent). However there is no need to investigate detailed 
transport through this complex mixed phase-space struc- 
ture, as the correlations give us a generic and quantitative 
handle on the energy diffusion with time. In this inter- 



mediate regime, dominated by the Poisson correlations, 
atoms prepared outside the trapping regions rapidly dif- 
fuse across the regions between them. Particles prepared 
in the trapping regions remain there. This regime occurs 
for t* ~ l/(Ke) 2 and corresponds approximately to the 
experimental results shown in FigEI 

Finally, at the longest timescales, there is the Cgi 
correction investigated in which is a long-ranged 

global-correlation family. Cgi results in an oscillation of 
the form — cospo^ and becomes dominant at the longest 
timescales. The oscillation is of the same period as the 
1-kick correlation, but is of opposite sign. This means 
that at the longest timescales, the minima in energy ab- 
sorption shown in FigUJa) become maxima in energy 
absorption; and vice-versa: the maxima become min- 
ima. Fig. 6(c) shows experiments tending towards this 
regime. The inverted peaks of the Poisson family are still 
in evidence, but a — cospo 6 oscillation is clearly super- 
posed. This is a somewhat counter-intuitive result since 
it implies that atoms initially prepared in the momen- 
tum trapping regions are the ones which at long times, 
for t* 3> l/(Ke) 2 , will absorb the most energy (there 
are no further reversals of this behavior at even longer 
times). 
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FIG. 7: Localization lengths of typical states for the double 
S-KP (e — 0.025, Ti = 0.5) corresponding to the three classical 
diffusion regimes investigated in the experiments in [Tll |. 

The classical analysis thus reveals that there are three 
distinct classical diffusive regimes occurring at three 
timescales. We can expect that the corresponding quan- 
tum behavior will depend on which regime is dominant 
when dynamical localization arrests the quantum mo- 
mentum diffusion. To investigate this, we now investi- 
gate how the form of the Floquet states varies in these 
different regimes. 

In Fig. 7(a) we show the localization lengths of the Flo- 
quet states for a weak kicking-strength. It can clearly be 
seen that the localization varies periodically as a func- 
tion of momentum, staying within the range 1 < L < 10 
for the majority of points, with the exception of a series 
of sharp cusp-like features at which the localization dra- 
matically falls. These location of these cusps exactly cor- 
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responds to the the "trapping momenta" , pe = (2n+l)n, 
predicted from classical arguments. The Floquet states 
centered in the trapping regions have widths of L ~ 0.01, 
much narrower than states localized on stable islands, 
which are also visible in this figure as regular strings of 
points at L ~ 1 . At the experimental values of K ~ 3 and 
e = 0.01 a similar behavior is produced, with the broad- 
est Floquet states having localization lengths of L ~ 60, 
while the narrowest have widths of L ~ 0.03, over one 
thousand times narrower. 

The Floquet states at the tips of the cusps have such 
low localization lengths that they are effectively pure 
plane-wave states (this can be further corroborated by 
evaluating the inverse participation ratio for these states, 
which indeed takes a value of almost unity) . It is thus un- 
surprising that the presence of these states corresponds 
to the classical trapping effect, as a quantum system pre- 
pared in such a state will have a vanishingly small over- 
lap with any other state and so will remain frozen (or 
trapped) in its initial state. It is important to note, how- 
ever, that this quantum trapping effect depends critically 
on the order of the two kick- periods - that is, whether 
the system is driven with a short-long kick-sequence or 
the inverse long-short ordering. 

This may appear surprising at first, since the Floquet 
states are periodic, with the same period T tot as the driv- 
ing, and this period is not altered by interchanging the 
order of the kicks. Although it is frequently neglected, 
however, it is important to recall that the Floquet states 
do have an explicit time-dependence within each period, 
and this is able to produce substantially different behav- 
ior 01 when the phase of the driving field is altered. To 
illustrate this, we show in FigOthe time-evolution of one 
of the localized Floquet states, which experiences £-kicks 
at times t = T\ = 1.90 and t = T to t = 2. As can be 
seen, the state has only a trivial time-evolution during 
the first time- interval (0 < t < Ti), since it is almost a 
plane wave and is thus approximately an eigenstate of 
the free Hamiltonian. The first kick at T\ causes the 
wave-packet to spread considerably in momentum space, 
before the second kick restores this broadened state to its 
original narrow form. Thus in this brief window of time 
between the two kicks, even the most localized Floquet 
states have a considerable spread in momentum. As a 
consequence, if the phase of the kicking field is shifted so 
that the system experiences the short-long kick-sequence, 
none of the Floquet states are sharply localized in mo- 
mentum at t = 0. We show in FigOHthe time-evolution 
of the system's kinetic energy when it is prepared in a 
momentum eigenstate in a trapping region. For the long- 
short kick-sequence this state projects onto essentially a 
single Floquet state at t = 0, and so its time evolution 
is trivial and its energy remains constant. For the case 
of the short-long kick-sequence, however, the initial state 
projects onto a number of Floquet states (Ea l3.2[) . giv- 
ing rise to a complicated quasi-periodic behavior arising 
from beating between the different quasienergies. 

In Fig. 7(b), we see the effect of increasing the kick- 




FIG. 8: Time evolution of a localized Floquet state in momen- 
tum space, N(p,t), for physical parameters: K — 2, e = 0.1 
and h = 0.5. Initially the Floquet state is sharply peaked at 
p = — 107T, in the center of a trapping region. The first kick 
at t — 1.9 causes the state to spread across a much broader 
range of momentum, until the second kick at T — 2 restores 
the localized state. 
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FIG. 9: Time-evolution of the energy of the double S-KP, 
for the same physical parameters as in Fig. 7a. The solid line 
shows the evolution of the system under the long-short kicking 
sequence, and shows little features. The dotted line shows the 
result of the short-long sequence, and exhibits a complicated 
quasi-periodic behavior. 



strength. In this regime we find there is an almost con- 
stant localization length for momenta in between the 
trapping regions, which are again signalled by sharp cusp- 
like structures. This indicates that the Floquet states 
are confined between the classical broken phase barriers 
in the trapping regions. Early studies indicate that the 
level statistics of the corresponding quasienergies are not 
pure Poisson in this regime, as would be the case for the 
standard QKP. 

In Fig 7(c), we see an inversion of the broad momentum 
modulation in Fig 7(a), similar to the reversal seen in the 
experiment. In this regime, the eigenstates localized in 
the trapping regions near p ~ (2n + l)7r/e are typically 
broader than those localized in between. 
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V. CONCLUSIONS 

We have presented a study of the Floquet states of 5- 
kicked particles pulsed with unequal periods, and used 
them to analyze experimental data on these systems. 
We conclude that the chaotic ratchet effect proposed in 
and observed experimentally in [t| is associated with 



asymmetric Floquet states localized around p — 0. We 
conclude also that the behavior of the localization lengths 
of the Floquet states for the double 5-kicked rotor broadly 
accompany the change over between the three distinct 
classical diffusion regimes investigated experimentally in 

m 



[1] G. Casati, B.V. Chirikov, F.M. Izraelev F.M., and J. Ford 
in "Lecture notes in Physics", Springer, Berlin 93 , 334 
(1979); B.V. Chirikov, Phys. Rep. 52, 263 (1979). 

[2] S. Fishman, D.R. Grempel, and R.E. Prange, Phys. Rev. 
Lett. 49, 509 (1982). 

[3] FL. Moore, J. C. Robinson, C.F. Bharucha, Bala Sun- 
daram, and M.G. Raizen Phys. Rev. Lett. 75, 4598 
(1995). 

[4] D. A. Steck, W. H. Oskay, M. G. Raizen, Phys. Rev. Lett. 
88, 120406 (2002) 

[5] W. K. Hensinger , H. Haffer, A. Browaeys, et al. "Dy- 
namical tunnelling of ultracold atoms" Nature 412, 6842 
(2001) 

[6] F. M. Izraelev, Phys. Rep. 196, 299 (1990) 

[7] T.S. Monteiro, P.A. Dando, N.A.C. Hutchings and 

M.R. Isherwood, Phys. Rev. Lett 89, 194102 (2002). 
[8] T. Jonckheere, M.R. Isherwood and T.S. Monteiro, Phys. 

Rev. Lett. 91, 253003 (2003). 
[9] P.H. Jones, M. Goonasekera, H.E. Saunders-Singer and 



D. Meacher, preprint arXiv:quant-phys/0309149 
[10] P.H.Jones et al, to be published. 

[11] P.H. Jones, M. Stocklin, G. Hur, and T.S. Mon- 
teiro, Phys. Rev. Lett. in press (preprint 
arXiv:physics/0405046 I 

[12] T. Dittrich, R. Ketzmerick, M.-F. Otto, and H. Schanz, 
Ann. Phys. (Leipzig) 9,1 (2000); H. Schanz, M.-F. Otto, 
R. Ketzmerick , and T. Dittrich, Phys. Rev. Lett. 87, 
070601 (2001). 

[13] S. Flach, O. Yevtushenko, Y. Zolotaryuk, Phys. Rev. 

Lett. 84, 2358 (2000). 
[14] P. Reimann, Phys.Rep. 361,57 (2002). 
[15] DL. Shepelyansky Phys. Rev. Lett. 56, 577 (1986). 
[16] B.G. Klappauf, W.H. Oskay, D.A. Steck, and 

M.G. Raizen, Phys. Rev. Lett. 81, 1203 (1998). 
[17] J. Ringot, P. Szriftgiser, J.C. Garreau, and D. Delande 

Phys. Rev. Lett. 85, 2741 (2000). 
[18] C.E. Creffield, Europhys. Lett. 66, 631 (2004). 



